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Abstract. We analyze the stress-energy tensor for the self-coupled Maxwell-Dirac system in the 
bilinear current formalism, using two independent approaches. The first method used is that 
attributed to Belinfante: starting from the spinor form of the action, the well-known canoni¬ 
cal stress-energy tensor is augmented, by extending the Noether symmetry current to include 
contributions from the Lorentz group, to a manifestly symmetric form. This form admits a tran¬ 
scription to bilinear current form. The second method used is the variational derivation based 
on the covariant coupling to general relativity. The starting point here at the outset is the tran¬ 
scription of the action using, as independent held variables, both the bilinear currents, together 
with a gauge invariant vector held (a proxy for the electromagnetic vector potential). A central 
feature of the two constructions is that they both involve the mapping of the Dirac contribution 
to the stress-energy from the spinor helds to the equivalent set of bilinear tensor currents, through 
the use of appropriate Fierz identities. Although this mapping is done at quite different stages, 
nonetheless we hnd that the two forms of the bilinear stress-energy tensor agree. Finally, as an 
application, we consider the reduction of the obtained stress-energy tensor in bilinear form, under 
the assumption of spherical symmetry. 


1 Introduction 

In classical field theory, it is well known that a conserved energy-momentum tensor for a physical 
system can be constructed by the Noether method, exploiting invariance of the Lagrangian 
density under space-time translations. However, the derived symmetry current using this method, 
the so-called “canonical” form uni, has the unfortunate drawbacks of being neither symmetric 
in its two tensor indices, nor gauge invariant. 

Many attempts have been made to rectify these problems, two of the most prominent being 
the Belinfante form [T], and the variational form from general relativity [22]. The basis of the 
Belinfante approach is to extend the invariance of the Lagrangian to include contributions from 
Lorentz transformations, so that the Noether symmetry current becomes that associated with the 
full Poincare group. In this way, the canonical term is symmetrized [23] . and extra “correction” 
terms are present, which are attributed to the spin contribution to the stress-energy. 

The variational approach uses the action principle to relate the variation of the Hilbert action 
of space-time to that of matter, then invoking Einstein’s equations to identify the matter part to 
the stress-energy tensor. The result is an expression for the stress-energy which is proportional 
to the functional derivative of the action of matter, with respect to the inverse metric. The 
presence of the metric ensures that the stress-energy tensor is manifestly symmetric. 

Regardless of the independent nature of their derivations, Goedecke pointed out [5] that in 
the limit of flat space-time, the Belinfante and variational forms of the stress-energy tensor must 
agree. The equivalence in the integral spin field case was proven by Rosenfeld m , and Goedecke 
provided evidence for equivalence in the half-integral spin field case via a series of examples, but 
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was not able to provide a general proof. Such a proof for the half-integral case was published 
shortly afterwards by Lord using the vierbein formalism. 

Much effort has been made to “improve” the stress-energy tensor, either by generalizing it 
beyond the Belinfante/variational forms, or by altering it so that it is compatible with a given 
theory. An example of the generalization aspect is the work done by Gotay and Marsden, who 
model the stress-energy in terms of fluxes of the multimomentum map across spacetime hypersur¬ 
faces [9]. The Gotay-Marsden stress-energy tensor naturally includes the spin “correction terms” 
present in the Belinfante formula, but in a more generalized fashion, as well as coinciding with 
the variational form in the presence of a space-time metric. Work has also been done by Callan, 
Coleman and Jackiw on making the cut-off dependent symmetric stress-energy tensor compatible 
with renormalized perturbation theory, by constructing appropriate counter terms in order to 
make it finite at arbitrarily large cut-off values [3]. The renormalization compatible stress-energy 
tensor is also compatible with an altered, but phenomenologically consistent, version of general 
relativity, and simplifies the currents associated with scale and conformal transformations in 
which the stress-energy appears. 

The work we present in this paper can be viewed as analogous to one of Goedecke’s examples, 
namely the coupled Maxwell-Dirac fields, but in an alternate formalism where the spinor fields 
are mapped bilinearly to a set of tensor fields. 

The background to this approach is as follows. Interrelationships between these “bilinears” 
can be derived [7] by taking their quadratic products, and expanding out the xV 7 spinor matrices 
in terms of the Dirac-Clifford algebra basis of 4 x 4 matrices. For the X — ^ case, there are sixteen 
real, gauge invariant bilinears, equal to the number of basis elements. Comparing this number 
with the eight real components of the Dirac spinors, and taking the gauge invariance into account, 
implies that there must be nine independent algebraic equations m- Furthermore, the bilinear 
map is invertible according to the spinor reconstruction theorem by Crawford j4j , so given a state 
described by the set of sixteen bilinear tensor fields, the Dirac spinor field is determined up to a 
phase. This bilinear reformulation is generalizable to the non-Abelian isospinor representations, 
at least in the SU( 2) doublet case, where an analogous set of Fierz identities involving sixty four 
bilinears can be shown to exist HU- 

There are many other identities involving the sixteen bilinears beyond the fundamental set, 
which can be derived by using the Fierz expansion along with known identities. The extension 
of this set of sixteen to other classes of bilinears, such as gauge dependent objects where some 
of the spinors are charge conjugated (i.e., % = ip c ), or contain derivatives, is a necessary part of 
the mathematical framework for the description of the Maxwell-Dirac equations in bilinear form 
m- A serious attempt at generalizing the fundamental set of nine Fierz identities to include a 
larger set of bilinear “currents” was presented by Takahashi [2T], although bilinears containing 
spinor derivatives were only discussed briefly. 

The bilinear representation of electrodynamics is compatible with the Maxwell-Dirac equa¬ 
tions, which model the self-consistent behaviour of an electrically charged fermion field interacting 
with its own electromagnetic field. The derivation of this system relies on the gauge covariant 
Dirac equation being invertible, with the gauge field, or vector potential A M , determined by 
the state described by the spinors, or equivalently, the bilinears. The inversion in the t/(l) 
(electromagnetic) case was originally considered by Eliezer [5], however he pointed out that the 
matrix required for the inversion had vanishing determinant, which discounted the possibility for 
a unique solution. More recently it was shown [2 that if the physical requirement that the vector 
potential is real is imposed, then the Dirac equation is indeed invertible, provided that a set of 
consistency conditions are satisfied. Incidentally, a formal inversion of the Dirac equation for the 
SU(2) Yang-Mills gauge field has also been achieved [IT] which extends to U{ 2) = SU( 2) x U( 1), 
but whether or not the inversion can be extended even further to the more physically relevant 
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E/(l) x SU( 2) x SU( 3) case is an open question. The set of equations obtained when the inverted 
Dirac equation is substituted into the Maxwell equations, along with the appropriate consistency 
conditions and Fierz identities, constitutes the Maxwell-Dirac equations. We shall not discuss 
the details of this system further here, other than to state the aim of this paper, which is to 
derive a bilinear form of the stress-energy tensor for this system, and explore how it reduces 
in the presence of symmetries. For further information on the bilinear current, or “relativistic 
hydrodynamical” formulation of the Maxwell-Dirac system, we direct the reader to [5D], and j!2j 
where symmetry reductions are also studied. 

The first topic we undertake in this paper is the derivation of the bilinear form of the Maxwell- 
Dirac-Belinfante tensor in Section [2] Following a brief derivation of the Belinfante tensor for a 
free Dirac particle in the spinor representation, we introduce the core concepts of the bilinear 
mapping and Fierz expansions. These ideas are then applied to the derivative-dependent bilinear 
terms that appear naturally in the free Dirac-Belinfante tensor, and an appropriate Fierz identity 
is derived, allowing us to rewrite the explicitly spinor-dependent form of the tensor exclusively 
in terms of bilinears. Note that the Fierz transcription to bilinears occurs as the final stage 
of the calculation, once the Belinfante tensor in Dirac spinor form has already been derived. 
A more detailed version of this derivation is relegated to Q3] The known tensorial forms of the 
electromagnetic interaction and Maxwell field stress-energies are then added to the free Dirac 
contribution, resulting in a manifestly symmetric and gauge independent bilinear form of the 
Maxwell-Dirac-Belinfante tensor. 

Section [3] presents an independent derivation of the Maxwell-Dirac stress-energy tensor, which 
in this case uses the variational form known from general relativity. Beginning with the La- 
grangian density for an electromagnetically interacting Dirac particle, and initially ignoring the 
Maxwell held contribution since we are mainly interested in the behaviour of the bilinear Dirac 
contribution, we convert it to its analogous bilinear form, using a contracted form of the Fierz 
identity obtained in Section [2J This is in stark contrast to the calculational method in Section [2j 
where the Fierz transcription took place last. A brief review of how the variational stress-energy 
is obtained is then given. Then, using the general relativistically covariant form of the bilinear 
Dirac Lagrangian, the variational stress-energy is obtained, and is found to be in agreement with 
the Maxwell-Dirac-Belinfante tensor. 

Finally, in Section SJ the bilinear Maxwell-Dirac stress-energy tensor is subjected to the re¬ 
strictions imposed by an example symmetry subgroup of the Poincare group, namely the spher¬ 
ical symmetry group SO( 3). This extends the discussion of the symmetry reductions of the 
Maxwell-Dirac equations from [12] , essentially providing the mathematical framework required 
to calculate the mass-energy and momentum fluxes corresponding to solutions obtained under 
this symmetry and others, in principle. The discussion of such Maxwell-Dirac solutions and their 
physical properties, as well as the treatment of other select symmetry groups, are intended to be 
presented in follow-up works. 


2 Maxwell-Dirac stress-energy tensor via Belinfante 

2.1 Belinfante tensor for a free Dirac particle 

The Belinfante stress-energy tensor is the fully symmetric analogue of the well-known asymmetric 
“canonical” form, which for free Dirac particles if] 

= -i , (1) 

^■For more details on mathematical conventions, refer to m and ECU- 
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that satisfies the conservation condition 


d l = 0. 


( 2 ) 


In fact, T is the Noether symmetry current corresponding to imposing the invariance of the 
free Dirac Lagrangian 


& = \ [^(^VO - (^V’)7 m V’] - mipi/j (3) 

under the translation group. The asymmetry arises from the fact that only translations are 
considered in the derivation of ©> which neglects rotational contributions to the stress-energy 
[5]. Since d3j) is invariant under Lorentz transformations, we can use the formula for the Noether 
symmetry current divergence 


0JS? .. d& 


= 0 , 


(4) 


to obtain the Noether current directly. A vanishing, manifestly antisymmetric expression is 
obtained 


i d a 7 ct , a^}4’) + = 0. 


(5) 


This can be interpreted as an antisymmetric combination of a symmetric tensor [23] . which we 
call 0 M ". A tensor form which reduces to the left-hand side of © upon antisymmetrization is 


QA- = T n» + 1 da - MY, ■ 

O 


( 6 ) 


Using m to replace the left-most anti-commutator bilinear term, results in the manifestly sym¬ 
metric combination 

= \(f^ + T^) - itu <?+ MY, , (?) 

where the spin density [IS], defined in this case as 

= -\^{Y,^ v }^ +MY, =o, ( 8 ) 

vanishes when taking into account the identity in the Dirac algebra 

{ 7 ^,a^} = 2e^ 7 57p- (9) 

We therefore obtain the form of the Belinfante stress-energy tensor for a free Dirac particle 

+ f^) = -i \ - (d^)Y^] , (10) 

in agreement with Goedecke 0. Note that the Belinfante tensor is conserved 


<9 /J 0 ,i " = 0, 


( 11 ) 


and is equivalent to the Noether symmetry current of the Poincare group. 
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2.2 Dirac bilinears and Fierz identities 

The bilinear mapping from the spinor fields x and ip to the corresponding set of tensor fields is 
of the form xTriP, where T^ = {/,, 7 ^, 75 } is the set of R = 1,16 basis elements of 
the Dirac Clifford algebra. When x = ip, we obtain the set of manifestly U(l) gauge invariant 
fields 


ipT R tp = (12) 

which are all real except for w, which is pure imaginary, as is the additional tensor *s llL ' = 
iplsa^ip. Note that a common alternate convention is to use the real pseudoscalar field w = 
ipry^ip- Alternatively, setting % = ip c , we obtain for the T^ = 7 ^ element, the gauge dependent 
bilinear vectors 


+ i = ipc^ip, 

(13) 

= Re['0 c 7 ^ip] = —{ip c l IJ ‘ip + ipx^ip 0 ), 

(14) 

= Im[^ c 7 At ' 0 ] = —(ipl^ip 0 — ip Cr y fl ip). 

(15) 


The matrix product of two spinors can be expanded in the T# basis, using the Fierz expansion 

111 1 1 

ipX= = ^(xV’)-f+-j;(X7MV’)7 /i + g(xo'^V’)o' /il/ -^(X757MV’)757 /i + -j;(X75V’)75, (16) 

R= 1 

where either \ or ip may contain partial derivatives. There exists a rich set of interrelationships 
between bilinear tensor fields, which can be derived by “breaking up” products of bilinears such 
as j^k v = ip'y ti (ipip)"/ 5 "f I/ ip into sums of different bilinears, by inserting (fl()l) . This Fierz expansion 
technique was used extensively in m, in the context of the reformulation of the self-coupled 
Maxwell-Dirac equations in terms of bilinears. 

2.3 Belinfante tensor in bilinear form 

Our current objective is to rewrite m in terms of Fierz bilinears, which means we need to derive 
a Fierz identity that expresses the spinorial object [ip') l ,(d ll ip) — {d^ipix^ip] in terms of bilinears. 
Therefore, we are led to search for Fierz expansions in which this term is likely to appear. One 
example is 

jui'tpidfj.ip) - (dfj,ip)ip] = - l -f{d^s va ) + ^{d^uj)k u - 

+ ^bP'lAd^ip) - (d^ip) 7„V>] - ^Suaiip^^id^ip) - (5 m ^)757 a ip] 

- 7^ a [ip') 5 (j v<T (di l ip) - {d^ip)^a u<T ip}. (17) 

There are at least three other bilinear products whose Fierz expansions produce the desired term, 
namely j^[ipl 5 {d^,ip) — {d^ip^hip}, k„(d^a) and k„(d^u ;). Their respective expanded forms, along 
with a much more detailed derivation of the Fierz identity, is given in appendix [B] Using these 
four identities, we can combine them to give 

[iplv^d^ip) - id^Xi,ip\ = (au)- 1 ^-^(d /u ,j a )(ojs I/a + a 7w) - k v [a{d^a) + w(d M w)] 
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( 18 ) 


+ \{8pk v ){a 2 +w 2 ) +j v {^{d^) - (d^)^] + - (d^ip)^]} 

+^(crs ua +^*s V a)[^l5l c ’{d^ip) - {d^)')5l' T ^] S j , 

which obviously still requires some more work. Using the Fierz identities derived in m 

[^{d^xp) - (d^ip] = ~(ct 2 - u 2 )~ 1 \j , '(d IJ ,k„)uj+ im v [d^n^a], (19) 

f^Ys^VO ~ (8pip) 75 ^] = -{p 1 - u 2 )~ 1 [f(dpk v )a + i m v {8^n v )oj], (20) 

we can replace the spinor terms on the second line of (fT51) with bilinears, but this still leaves the 
spinor terms on the third line. After a straightforward, but tedious, set of Fierz manipulations, 
we obtain the desired identity 

^(as„ a + u}*s ua )[tp-/ 5 -/ a {dpip) - (a /i V’)757 <T V’] = “<MV>7 Ad/ti’) - (d^huip] 

+ T >( cr2 - u} 2 )~ 1 [j a (d^k cr )(a 2 +w 2 ) + 2im CT (5 M n (T )crw] - |(ct 2 + w 2 )(5 /x fc J/ ) 

4 o 

+ ^k u [a(8pa) + w(<9 m w)] - Udp,j a )(a*s va +ws ra ) + \f[<j{8p*s va ) + w^s^)]. 
o o o 

( 21 ) 


Substituting into (HI and rearranging yields 

[ihfuidni)) - (S^b^] = (cw) -1 |-^( 0 ' 2 - w 2 ) _1 [j CT (5 M fc CT )(CT 2 +w 2 ) +2i m a (8p,n a )au}\ 
+ +u 2 )(8 ll k u ) - ^k v [a{8piT)+uj{d li w)\ + [<r(^*s w )+w(9 /f s w )] 

- ^(^j <T )(o’*s w + ws^)| . ( 22 ) 


This is entirely in terms of bilinears, but we would like to go further and eliminate the rank-2 
terms, s M „ and *s pv . Using the known Fierz identities for the replacement of these terms O 

= (a 2 - w 2 )-Ve^ CT - u8 pvpa )j p k G , (23) 

V" = (a 2 - w 2 )- 1 (we^ ff - aS^ pa )j p k a , (24) 

where we define the partially antisymmetric object 

8 pupa = _ r f" r j'P) i (25) 

the rank-2 dependent terms in (1221) become 

Y^j a W(d P *Sucr) + u{d^s va )] = ^(cr 2 - w 2 ) -1 {2crwfc y [w(5 M cr) - a(d p u)\ 

+ 2\aue va p t f(dpj p )k e + j„j a (8p,k a )(a 2 +w 2 )} - -^{8p,k v )(a 2 + u 2 ), (26) 

and 


5i 


5i 


- + us va ) = (<r -w ) j -^3vf{dp,k'j)((j + OJ ) -—aoje vape (dpj (7 )j p k' 


12 " 
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( 27 ) 


- + W 2 )[w(<9 m w) - cr(<9 p cr)] j , 

giving us the final form of our identity 

- (3^)7^] = ( a 2 - ui 2 )~ 1 {k^[u}(d IJ .a) - cr(<9 p w)] - ie vape (dpj a )j p k e 

- i j v m a {d^n G )}. (28) 

Substituting into m and relabelling some indices for convenience later on, we obtain the bilinear 
form of the Belinfante stress-energy tensor for a free Dirac particle 

O^i/.d = ^(ct 2 - w 2 )^ 1 {-i[fc / _ 1 (w9 l/ cr - ad v (jj) + k v (ujd p a - ad p u))] 

- j*k K [e u paK (d p ,j p ) + e tl paK {d v j p )\ - j p vF(d„n a ) - j v m a (d p n a )}. (29) 

The Fierz identity (1281) is of central importance to our current work on the stress-energy tensor, 
since it permits an explicit mapping from the Dirac spinor fields to the physically equivalent set 
of bilinears, in the form in which the spinors appear in the Belinfante tensor m- Of note is 
that none of the elements of the bilinear set CE2D which we have chosen to work with contain 
any internal derivative operators, so all gradient terms appear explicitly in (EH1) . Additionally, 
(1281) , along with Cl and d, are examples of a large family of related Fierz identities, and only 
represent the tip of the iceberg of such relations m- 


2.4 Maxwell-Dirac Belinfante tensor 

The full Maxwell-Dirac stress-energy tensor is 

©/ii^,MD — O“b © P iy,int “b (30) 

where ©^int and ©^em are the interaction and Maxwell field contributions respectively. The 
Maxwell contribution has the well known form 

© W m = - F pa F v a , (31) 

consistent with our metric signature (-1 -), and the interaction term is 

©^,int = -jii^Av F ji/Ap), (32) 

where the electromagnetic vector potential A p can be replaced by the gauge independent analogue 
B p using the definition from [T 2 ] 

Ap = Bp + -^(ct 2 - u 2 )~ l m a (dpn a ). (33) 

Zq 

The gauge dependent bilinear terms in (15^1) cancel out the corresponding terms in (EH1) exactly, 
so the full Maxwell-Dirac Belinfante stress-energy tensor is 

Qp V ,MT>=^((r 2 -oj 2 )~ 1 {-i[kp(ujd v a-ad u uj) + k v (ijdpa-adpU})]-j <T k K [. e„ p(TK (dpj p )+ep paK (d v j p )]} 

+ | {j»B v + j v Bp) + -^ v F ap F° p - F prr F v \ (34) 

which is manifestly symmetric and gauge independent. 
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3 Maxwell-Dirac stress-energy tensor via general relativity 

3.1 Bilinear form of Dirac Lagrangian 

We will now derive the bilinear form of the Maxwell-Dirac stress-energy tensor again, this time 
using a completely different method. The approach outlined in section [2] involved the use of Fierz 
identities to convert the spinorial form of the Belinfante stress-energy for a free Dirac particle 
m into bilinear form m, to which the known tensor forms of the interaction and Maxwell 
contributions were added, yielding (1551) . 

This time around, we convert the Lagrangian for an interacting Dirac particle 

% = - (<9^)7^] - fni/np - (35) 

into bilinear form, then use the definition of the stress-energy tensor from general relativity to 
obtain our result, which should in principle agree with (1511) . Note that this method from general 
relativity was used directly on the spinorial form of the Lagrangian ([3]) by Goedecke [3], who 
then demonstrated its equivalence with the Belinfante stress-energy tensor HDD- We are pursuing 
a similar equivalence demonstration, with our focus being on the bilinear formalism. However, 
in contrast with the calculation of the previous section where the Fierz bilinear transcription 
was performed at the end, once the spinorial Belinfante tensor m had been obtained, here 
we transcribe at the start , and perform the entire variational calculation with the bilinear field 
variables. The Fierz identity required to rewrite (1351) is obtained by simply substituting the 
contracted form of (1531) . 

[^(d^ip) - (<9^)7 / V] = (c 2 - w 2 ) _1 {/c A1 [w(a At cr) - cr(<9 M w)] - ie tUTpt (d IJ .j (r )jpk e 

- i j p 'm a (d p n< J )}. (36) 

Applying the definitions a = j p and using the definition of B p (1551) . we obtain 

= i(cr 2 - u 2 )~ 1 {ik p [u}(d p a) - cr(<9 p w)] + e paKT (<9 p j CT )j K fc r } - ma - qj p B p . (37) 


3.2 Variational form of the stress-energy tensor 

The total action for the gravitational field in the presence of matter is ;22j 


S = 


Sr 

167tG 


+ Sm, 


(38) 


where Sm is the action for matter fields (mass-energy). Sh is the Hilbert action, defined as 


S'h = / d 4 x 


(39) 


where g is the determinant of the metric g pui R is the Ricci scalar, and d 4 a; is the invariant volume 
element. The variation of the action with respect to an arbitrary tensor field 
the general form 


IJk u 1 ...u l takes 


SS = [ d 4 a^<5$ 
I <5$ 


(40) 
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with contraction over the indices implied. The term 6S/5$> is called the functional derivative of 
S with respect to the tensor field $. Of main interest in variational theory are tensors < f>o which 
extremize the action, so that SS = 0, and hence 


SS_ 

<5$ 


= 0 . 


(41) 


$0 


Extremizing the variation of the Hilbert action (1391) with respect to the inverse metric leads to 
Einstein’s equations in vacuum 


‘ = R „„ _ Igg = 0. 

yj^g Sg pi/ M 2 

Likewise, extremizing the gravitational action in the presence of matter 

js -iS +is “= 0 - 

and equating the corresponding functional derivatives, yields 


so that 


1 SS 


1 


\J~-g 8g lllJ 167 tG 


Rfw 2 R'9^1' 


1 SSm 

V~9 S g 


fllS 


= 0 . 


Comparing with the well-known form of Einstein’s equations in the presence of matter 


R/iv c\R9/iv — 8tv( 1 / /j,/. 


we can see that the energy momentum tensor is of the general form 

rr 2 SSm 

‘“'--yFgSgi"- 


(42) 


(43) 


(44) 


(45) 


(46) 


3.3 Variational Maxwell-Dirac stress-energy tensor 

The variation of the electromagnetically interacting Dirac matter action, ignoring the contribu¬ 
tion of the Maxwell field itself for now, is given by 

SS = J d A x (6y/=g& + y/=gS££) , (47) 

where the Lagrangian is given by (EH), but with appropriate modifications to make it covariant 
in curved space. Since the invariant volume element d 4 x and yj—g are scalar densities of weight 
— 1 and +1 respectively, we must arrange for the Lagrangian to be manifestly a scalar. Notice 
that E3 contains a term dependent on the Levi-Civita symbol with upstairs indices, which is 
of weight —1. This implies that we should make the replacement 

e paKT — =e paKT . (48) 

V^9 

In order to deal with the bilinear four-vectors we must introduce the vierbein fields [Mi, which 
locally relate the curved metric to the flat one 

9/iv — &/1 6 V 7]abi (49) 
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where Greek and Latin indices label curved and flat spacetime components respectively. The 
gamma matrices are modified such that 


TV — l a ' {la, lb] 2?y a f,, 

so the bilinears are now 

j/x = e^ja, ja = Iplalp, 

*V = e M a /c a ; k a = ipiniaip- 

The variation of the square root of the negative metric determinant is m 


= - 2 ^~99^6g tlv . 


(50) 

(51) 

(52) 

(53) 


Noting that h = y/—g, where h is the vierbein determinant, we can use the variation of (147)1) to 
alternatively write this as 

Sh = -he p a (6e p a ), (54) 

implying the reciprocal variation 

Sih- 1 ) = h~ 1 e p a {8e p ‘ a ). (55) 

In curved space, the Levi-Civita term in (1371) becomes 

h~ 1 e p,TKT {dpe^jaj^r + h~ 1 e paKT e a a {d p j a )j li k T . (56) 

Introducing the covariant derivative causes the first term to vanish, due to the tetrad postulate 


S7 p e„ a = 0 . 


Expanding out all of the vierbein fields in the second term, we find that 
e paKT e a a {d p j a )j K k T = e abcd {d a j b )j c k d , 


(57) 


(58) 


which implies that for any curved coordinate components, this term is always equal to the 
flat spacetime version, so it is automatically covariant. We find that the covariant bilinear 
electromagnetically interacting Dirac matter Lagrangian has the form 

3? = ^{a 2 -Lo 2 )~ 1 [ig pcr e cr a ka(u}d p a—ad p Lo)+h~ 1 e paK ' T e cr a e K b e T c (d p ja)jbkc}-^ r ncr—qg p,T e cr a jaB p . (59) 
The variation with respect to deformation of the vierbein field is 

S3? = ^(a 2 - u, 2 )- 1 {i5(g prT e a a )k a (ujd p a - ad p u) + e p ™ T 8(h- l e a a e K b e T c ){d p j a )j b k c ] 

- q5(g p °e a a )j a B p . (60) 

From the variation of (147)1) . we find that 

% p<T e CT Q ) = 2 (<50 + (, 5e a b )e pb e a a . (61) 
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Using the fundamental vierbein property 


e p a e p b = 5^, (62) 

e^iSe^b) = -(Se^t, (63) 

we find that the first and last terms in m are 

i5(g pa e a a )k a (ujdp(j — ad p oj) = —ik p (ujd a a — ad a oj)(5e IJ ' a ), (64) 

- q5(g pa e a a )j a B p = qjpB a (5e p a ). (65) 


Following a similar process, we find that the second variational term is 
e paKT 5(h~ 1 e a a e K b e T c )(d p j a )j b k c 

= h- 1 [-ep. a e p ™ T e rJ d (dpj d )j K k T +e P p™(dpj a )j a k K +e p % K e a b (dpj b )j a k K +e p ™pe a b (dpj b )j K k a }(6e p a ). 

( 66 ) 

Gathering the deformed terms together, we can write the variation of the Lagrangian as 

&& = Q(<7 2 - - ad a u) + h- l [-e p a e p ™ T e a d {d p j d )j K k T 

+e p fl ™(d p j a )j„k K + e p \ K e a b (d p j b )fk K + e p ™ p e a b (d p j b )j K k a }} + qj p B a ) (5e p a ), 

(67) 


with the associated action variation being 

5S U = Jd 4 x^g (-e p a ji(<r 2 - uj 2 )~ 1 [ik p (uid p a-ad p u})+e p,TKT (d p j C7 )j K k T ]-ma-qj p Bp | 
+ i(a 2 - cc 2 )- 1 !-! k p (cod a a - ad a u) + e a d (d p j d )j K k T + e p p ™(d p j a ) Ja K 


+ e p \ K e a b (d p j b )j a k K + e p ™ p e a b (d p j b )j K k a }} + qj p B a ) (6e p a ). (68) 

From the general form of the action variation (TO , a relationship between (1S51) and the stress- 
energy tensor can be obtained m 

5Sn = J d 4 x^gu\ a Se x a = ^ j d 4 Xy^gT in '5g pl/ , (69) 

which implies that 

u P a = ^(T pX e Xa + T Xp e Xa ). (70) 

Recognizing T pl , as symmetric gives 

T pv = + e vaU p a ). (71) 


Identifying the contents of the external parentheses in (1681) with u p °, we obtain for the stress 
energy tensor 

T pl , = - ripv j^(cr 2 - w 2 ) -1 [i k p (ud p o - ad p io) + e pcrKT (d p j a )j K k T \ - ma - j 
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+ j(a 2 - uj 2 ) l {-i[kp{ud v a - ad u u}) + k u {ud p a - ad p u:)} - 2r] flu e paKT (d p j a )j K k T 
+ z P n crK ‘{dpjv)jcrk K + e pa P K (d p ,j (T ). j v k K + e paK p(d p j <T )j K k v + e p v aK (d p j p )j a k K , 


+ e pa v K (d p j a )j lt k K + e p °\(d pJa )j K k p } + \{j p B v + j u Bp), (72) 

where we have evaluated at flat spacetime. This is manifestly symmetric, but it requires some 
additional manipulation before it more closely resembles the Belinfante form (|34|) . Consider the 
U(l) gauge covariant Dirac equation and its Dirac conjugate 

iy cr (<9 CT V’) “ 97 <T ArV’ - mip = 0, (73) 

i{d a ip)'y C! + qil>'y a A G + mip = 0. (74) 

Left and right multiplying these equations by ip and ip respectively, then subtracting the second 
from the first and rearranging, gives 

m.a = ^ ^^{daip) - ( d a ip)"i a ip\ - qj a A a (75) 

Applying the Fierz identity (E!6l) and the B p definition (1331) . this becomes 

-Trier = - i(a 2 - cc 2 ) -1 [ik p (ujd p a - crd p u) + e paKT {d p j a )j K k T \ - qj p B p , (76) 


causing the rj pl , dependent term in (1721) to vanish. Now consider the combinatorial identitj@ 

- e p paK -{d„j p )j a k K - £v paK (d fl j p )j a k K 

= —2rj p , u € paKT (dpj^jnkr + e p p , aK (d p j u )j a k K + e p % K (<9 p j a )juk K + e paK n{d p j a )j K k v 
+ c P i/ ,TK '(d p j tl )j CT k K + e pcr v K {d p j (7 )j l j,k Kj + e paK i/(d p j cr )j K kp, (77) 

which can be used to obtain the final form of the variational stress-energy tensor for Dirac matter 

Tp V p = ^{a 2 -u 2 y i {-\[k p {ud 1J (T-<jd v u)+k v (ujd p ,(j-<jd IJ ,u)\-j„k K [ev P ' TK {d p ,jp)+e p , p ' TK {d v j p j]} 
+ | (jpBv + j„Bp). (78) 

Comparing with the Belinfante tensor (1341) . we find that they agree 

Tpv, MD = 0/xi/,MD, (79) 

when the gauge field stress-energy m is included on the left-hand side. 

4 Symmetry reduction of the Maxwell-Dirac stress-energy 
tensor 

The Maxwell-Dirac equations in the bilinear formalism are in general, a very complicated set 
of self-coupled partial differential equations. An application of the present construction of the 

2 This follows from the 5 term cyclic identity V a e paKT + V T e apr7K + ■ ■ ■ = 0 which holds for the Levi-Civita 
tensor multiplied by any contravariant vector quantity. With the role of V a played by the Kronecker rill (for fixed 

/3), this yields 1771 1 after contracting with ri pa S^ (dpj^j «At , and rearranging indices appropriately. 
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physical stress-energy tensor of the system in terms of bilinears, is therefore to provide a rep¬ 
resentation of the conserved rest mass of possible solutions (via the spatial integral of Too for 
example). 

For any meaningful solutions to be derived, it is natural to consider reduced forms under 
the imposition of special symmetries. The reduction of the bilinear form of the Maxwell-Dirac 
system under several examples of subgroups of the Poincare group was discussed in [12] . We 
therefore choose one of the most important of these subgroups to work with here, namely SO{ 3), 
which corresponds to spherical symmetry. In particular, we shall demonstrate how the bilinear 
form of the Maxwell-Dirac stress-energy tensor reduces, given the restrictions imposed by this 
subgroup. The treatment of other symmetry reductions, such as cylindrical symmetry and the 
A 3,10 subgroup from {IS], shall be left for future work. 

Under spherical symmetry, scalar fields (a, w,etc.) have the generic form 


</> = f(t, r) 

and vector fields (j^, fc M , etc.) have the form 

/ \ 

= xg{t, r) 
yg(t,r) 

\ zg(t,r) ) 

where the invariant 


(80) 


(81) 


= y/x 2 


y* 


(82) 


is simply the spatial radius. We showed in m that the spherically symmetric forms of our 
bilinear vector and axial vector fields are 


f = 


( Ja ) 


/ 

rjb 

\ 

xjh 



(■ x/r)j a 


yjb 

, AC — 


(' y/r)ja 


\ z ib y 


l 

( z/r)j a 

/ 


B» = 


( B a \ 

xB b 
yB b 
V zB b ) 


(83) 


where the vector potential functions are 

B n = ±^(oyw - auj r ) - maj a [q{a 2 - uj 2 )] 


B b = 


T —(cr t w -(not)- maj b 
2 r 


[g{o- 2 -u 2 )\ 


-1 


(84) 

(85) 


Here, we are using a condensed derivative notation dtc r = at, and so on. Note that the effect of 
the symmetry reduction has in this case, reduced the components of the four-vectors to the set 
of coefficient functions j a , j b , er and w, which are themselves further constrained by higher-order 
nonlinear PDEs in the Maxwell-Dirac system. The coefficient functions k a and k b have been 
eliminated through the use of the Fierz identities 


a u 


j^ = 0 . 


( 86 ) 

(87) 


It is straightforward to show that the Levi-Civita terms in the stress-energy vanish in this sym¬ 
metry case 


Me u ^ K (dJ p ) + e^ K {d v j p )} = 0 . 


( 88 ) 
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The form of the stress-energy tensor we are dealing with is therefore 


V,md = - ^(cr 2 - w 2 ) 1 [k li (ujd v a - ad v bj) + k v {ujd^a - ad^ix)} + | {jfj,B v +j„B fl ) 

+ (89) 

For the SO( 3) symmetry, the components of the electromagnetic field strength tensor are 

F 0i = XiF a (t,r), (90) 

Fij = e ijk x k F b (t, r ), (91) 


where the Maxwell coefficient functions are 

F a (t,r ) = (1 /qr){cr 2 - w 2 ) _2 {-2 m[aj a (aa r - ww r ) + raj b (aa t - uu> t )] ± i[crw(cr 2 - of + w 2 - of) 
+ (cr 2 + w 2 )(o t w t - <r r w r )]} + (l/gr)(o 2 - u> 2 )~ l [m(a r j a + aj a , r + ra t j b + raj bjt ) 

± (i/2)(cr tt w - autt - cr rr ui + auj rr )\, (92) 


representing the electric field form, and 

Fb{r) = ± 2^' (93) 

representing the magnetic field form, which happens to be that of a monopole. Treating the 
H = is = 0, H = 0, v = i and /i = i, v = j cases separately, we find that the respective 
components of ([89]) are 


Tr 


00, MD 


— T a + T , 


Xi 

rji t rji 

-^Oi.MD — l-bi 

r 

XjXn „ _ 

Tij, md = + SijF , 

where the energy density of the Maxwell field is 

r\Fl + F^ 

2 

and the other functions are defined as 


T a {t, r) = (a - w 2 ) \ ±- [j a (u r uj - au> r ) - rj b {a t w - aw t )\ - mcrj 


2 

a ( 5 


T b (t, r) = (cr - w 2 ) ±- \j a (a t u] - aui t ) - rj b (a r u) - au; r )] - mcrj a rj b 


T c (t,r) = (a -ui ) ±- [j a (cr r uj - au r ) - rj b (a t u] - au> t )} - mar j b - 2 T. 


(94) 

(95) 

(96) 


(97) 

(98) 

(99) 

( 100 ) 


5 Conclusions 

The central aim of our current work was to obtain a form of the stress-energy tensor for the 
self-coupled Maxwell-Dirac system, in terms of bilinears. The motivation for working in the 
bilinear formalism in general, is given by the fact that spinor fields do not correspond directly to 
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observables, due to their dependence on an arbitrary unobservable phase factor. On the other 
hand, the system of Dirac bilinears corresponds to “density fields”, or “probability distributions”, 
which are observable according to a fundamental postulate of quantum mechanics. However, the 
attractiveness of the new physical insights that could be gained by studying the fully non-linear 
quantum electrodynamics, in terms of the observable tensor fields given by bilinears [50], is 
offset somewhat by the sheer complexity of the self-coupled system. Reducing this system under 
a symmetry subgroup of the Poincare group can simplify it to a point where it is mathematically 
tractable, and solutions can potentially be obtained. Such solutions would have associated mass- 
energy and momentum fluxes, which can be calculated directly using a bilinear-dependent form 
of the symmetry reduced stress-energy tensor. 

In sections [5] and [3] we demonstrated that bilinear forms for the stress-energy tensor (0^ 
and Tj lv respectively) can indeed be calculated, by applying Fierz identities to the spinor terms 
appearing in the Belinfante and variational general relativistic calculational schemes. In the 
Belinfante case, the Fierz mapping was applied to the spinorial Belinfante tensor, and in the 
variational case, it was applied to the spinorial Lagrangian prior to the vierbein deformation. 
Despite the fact that these two methods are quite independent, with very different starting points 
from the point of view of the Fierz bilinear transcription, they are in agreement in accordance 
with Goedecke’s conjecture [5] and Lord’s subsequent equivalence proof mi- This mathematical 
consistency adds weight to the validity of the bilinear representation of the Maxwell-Dirac system, 
and our corresponding bilinear stress-energy result. 

However, there is a point of view from which this automatic agreement is somewhat surprising. 
When taking into consideration the details of the functional Jacobian corresponding to the spinor 
to bilinear mapping, one would expect there to be extra constraint terms entering into the 
bilinearized Lagrangian, with the lack of such terms in m leading to a disagreement with the 
Belinfante tensor in the bilinear representation. A transcription of spinor electrodynamics into 
gauge invariant quantities in this spirit, has been given in the functional formalism by Rudolph 
and Kijowski HE ns]- In their bosonic transcription, Green’s functions are given as functional 
integrals in whose integrands there are always additional accompanying field-dependent factors, 
and so an effective bosonic, local, purely Lagrangian formulation is not obtained. The details of 
the agreement between 0 M1/ and for the bilinear case, although highly encouraging, remains 
a matter deserving of further study. 

Putting these technical concerns aside, we then turned to an example to demonstrate how 
the bilinear stress-energy tensor is reduced under spherical symmetry, using the generic SO( 3) 
invariant forms for scalar and four-vector fields discussed in m ■ We found that the stress- 
energy components could be described in terms of three functions (fMl) - (ll()()l) corresponding to 
the interacting Dirac matter contribution, as well as a single function EB, corresponding to the 
energy density of the Maxwell field. 

The bilinear form of the Maxwell-Dirac stress-energy tensor (1351) may be easily applied to 
solutions of the symmetry reduced Maxwell-Dirac equations, such as considered in m , by per¬ 
forming the corresponding symmetry reduction on T Ml/ , as done in section[5]for the SO( 3) group, 
and directly substituting the solution fields. We intend to investigate both numerical and closed- 
form solutions explicitly, for the static spherically symmetric and the algebraic splitting groups 
-Pi 3 ,io |35j, [13, in follow-up works. 
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A Dirac identities 


{ 7 ^, 7 "} = 2 »T 
[ 7 ^, 7 "] =~ 2 i 

7 5 = 75 = -(i/4!)e M1 ,p C r7 M 7 I/ 7 P 7 IT = i7°7 1 7 2 7 3 = -170717273 
75=4 

{ 75 ,7 P } = 0 

[ 75 ,^] =0 
7 ^ 7 ^ = ^ - ‘\(7 pl ' 

7 p 7m = 4 
7 p 757m = 4 75 

7 M 7 - 7 A = ?? ^ 7 A + ^Ayx _ ^A^ _ ie ^Aa 757CT 

= -27 m 


= 2757 ^ 

'f'f'f'f = 77 ^'"r 7 CTe + - 77^77^ 


vrf 






irf e a‘ rv 


+ i^V^ + i^tr"'* - ie pi/CT£ 7 5 
'fa*"' = - \rf v ^ + e^ £CT 757 CT 

= ir? - 7 M _ ir] ^Y + e /— 757ct 


7 m ct < 7£ 7 1/ = i77 £ V CT 




77 £ly cr^ CT - Tj rTL '(j^‘ e - e° 


L-, 1 • crevX-, _n 

75 + ie 75^ a 


7^7* = 0 
^"7m = -3i 7 ", 
a fJ ‘ v, y p ') li = 3irfP + tr iyp , 

= r\ vPr f - rTY + 

7 p aV/* = — 3i7„, 

7 P 7 p 0 V At = 3i(5„ p - cr I , p , 

7 p cr pT (7^ /i = $„ T 7 P - <W + i7-e Kpr ' T 757a, 

- e ApCT£ e A pi/T = r] pp rf v rf T - r) Pil rf v rf T + r] pv rf T rf p - 7 f v rf T rf p + r] pT rf p rf v 


B Derivation of the Belinfante Fierz identity 

Here we supplement section 12.31 with a more detailed version of the derivation of (1281) . 
Fierz expansions containing the term we want to solve for, [(d fJ ,ip)"/v' l P ~ V’7i/(A' ! /’)] are 

jdV’tAVO - (^V’)V’] = ^(9 M j !T )w - ^j CT (A s ^) + 

+ ^tr[V’7«/(^V’) - (AV’b^V'] - ^•w[b757 <T (d/^) - (3 M ^)757 a ^] 

- -^[-tp^avaid^ip) - (AVO750W0], 

bhbbAVO - ( 9 ^) 75 ^] = ^(Ai'T'W - ^T(V<w) + ^OV)^ - 


( 101 ) 

( 102 ) 

(103) 

(104) 

(105) 

(106) 

(107) 

(108) 

(109) 

( 110 ) 
( 111 ) 
( 112 ) 

(113) 

(114) 

(115) 

(116) 

(117) 

(118) 

(119) 

( 120 ) 
( 121 ) 
( 122 ) 

(123) 

(124) 


The four 


(125) 
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+ ^bPlA 9 ^) - vV’] - 757 ct (5 m V’) - (3 / xV’)757 a V’] 

- ^k <T [ipa V(T (d tl il>) - (d^owV’]) (126) 

ku(d^a) = ^(d^K) - ^{f*s VC7 ) + [^75 (d^) - (0^)75^] 

+ ^s j/(T [V’757 ct (^V’) - (5 m V’)757' t V'] - ^[^^(<9^) - (<9^V0°W^] 

- ^[HA 9 ^) - (<9 J uV’)7i'V’]> (127) 

kvid^w) = ^(dpK) - ^d /J ,(j a s„ cr ) + ^j„[ip(d M ip) - (d^ip] 

+ ^*Sva[iH5l a (d^) - {d^)lsl a ^] - ^k a [ip'y 5 a va {d, jl ip) - (d^bscwV'] 

- \aVI>lv{d^) - {d^huip\- (128) 

Combining these equations gives 

Vl>lv{d l xi>) ~ (0^)7^] = (aw)" 1 (^-^{d, Jl j cr )(ujs ua + cr*s ua ) - k v [a{d^a) + w(0 M w)] 

+ \i 9 nk v ){(J 2 + w 2 ) + j„{w[V>(0 M ^) - (0 M V#] + o#75(0/WO - (<9^V ; )75'0]} 

+^(as yCT +w7w)[^757 CT (<9^) - (0^)757^]^ , (129) 

which obviously requires more Fierz manipulation, since there are still spinor terms present. 
Using the Dirac identities (U18I) - (I120[1 . we obtain the additional Fierz expansions 

s va [ip'y 5 'Y a (d^ip) - (<9^)757 a ^} 

= | l -a(d li k v ) - + ^j a {d,j*s va ) - jjr(0 f J‘ T )*s I , o - + ^>[^75(9/iV’) - (3^)75^] 

- - (d^CTvcrlp} + y LJ[ijry u (0^) - (d^ivip], (130) 

*SvabP r Y5l cr (d tl ip) - (0^)757^] 

= ^{dpkv) - ^{dnw)k v + ij CT (< dy.s va ) - i( d,j,f)s VC r + ^jvbK 9 ^) ~ (0^)^] 

- jU^hh'sow^V’) - (0^)7 %a va ip] + ^a[^7^(0 M VO ~ (0^)7*# (131) 

Combining these expansions into the form they appear in (11291) . we get 

^(crs„ a + u}*s„ a )[4>'y 5 'y a (d^) - (d /i V’)757' T '^] = ^M^^a) + w(^w)] 

- Y^(a 2 +uj 2 ){d^K) - ■^{ 9 i_ L j a )(a*s vcr + ujs„ a ) + ^f[u{d^s va ) +u(5 fl s w )] 

+ ^>(a 2 - .r^(W(a 2 + w 2 ) + 2im <T (0 M n cr )<Tw] 

- ^-k a {ablxjyci.d^ip) - (<9^)07^] + u#750w(0 M VO - 
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3 — — 

- -CTUj[^ v {d^) - (< 9 ^) 7 ^], 


(132) 


which itself contains terms requiring further Fierz analysis. Using the Dirac identities (11211) - 
(11231) . we find that the expansion of these terms is 

- (d^VOcw^] = ^(d^j a )*s va - ^j a (9^*s va ) + ^{d^a)k v - ^<j{d, Jt k v ) 

+ jjvbhs (^VO ~ (<^V075^] - jr^o# 757 <T (5 / xV’) - 


jU>[ip'Yvidnip) - (d^ip) 


(133) 


k a [iri^v<r{d^il>) - (d^VO75<W^] = i(9 M j' T )s l/(T - + ^{d^uj)k u - yu^*;,,) 


y^W^) - (^V’)V’] - ^*s^a[^757 CT (^^) - ('9 m V’) 757 <t ^] 


3i 




^ “ ir ivx-'nH'J K'-'l-i 1 , 

Again, combining these terms into the form in which they appear in (11321) gives 


(134) 


- — k a {a[ipa va {d il ii) - (d^a^ip] + ui[ip'YBO’va(d tl ip) - (0^)75owV’]} 

= ^(crs^ +w*s 1 /ct )[V>757 <t ( 5 m ^) - (^V’)757 <t ^] + ~ (<^V 0 lA] 

- ^ji/(o- 2 - u; 2 )~ 1 [j' 7 (d f _ l k a ){a 2 + u 2 ) + 2im a (d^n a )au)\ - ^(cr 2 + u 2 )(d tl k v ) 

+ ^k v [a{d^a) + w(d M u;)] + - ^(d^Xo^cr +usj, 

(135) 

which when substituting into (11321) and rearranging, gives 

^(as yCT +wTw)[X757 CT (^X) - (<9 / /0)757 <t '*/’] = ~™bP / yv( d »' l P) ~ (<^X )l^\ 

+ tXX^ 2 - ^ 2 )~ 1 [j a (d^k cr )(a 2 + w 2 ) + 2im CT (<9 M n (T )cru;] - |(cr 2 + w 2 )(5 /x /c J/ ) 

4 o 

+ |fc J/ [cr(S /i c7) + widftU})] - \{d iJ f){a*s VCT +ws w ) + + w^s^)], 


a pure bilinear tensor expression. Now, using the Fierz identities derived in m 
[ip(d,j,il>) - {d^ip)^] = -{a 2 - u 2 )~ l [j"(d^k v )u + im"(d^n^al, 

[V’7s(^/i‘0) - (d^) 75X] = -(^ 2 - w 2 )~ 1 [T(d M M cr + im v (9 M n„)a;], 

and combining them into the form in which they appear in (11291) . we get 

w[V’(5 /i V) - (9^)^] + (t[4> 75 (^V’) - (^^)7 5 -0] 

= -(cr 2 - w 2 ) _1 [j‘ T (5 M fc cr )(cr 2 +w 2 ) + 2im CT (<9 AI n cr )cru;], 
which along with (1 1361) . can be substituted into (I129D to give 


(136) 

(137) 

(138) 


(139) 


[ir/Ad^tp) - {d^'y^] = (crw) <{ -- j v {a~ - w ) [j a {9,j,k a ){a + w ) + 2i m a {d, J ,n a )aoj\ 
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(140) 


+ + co 2 )(d p k u ) - j^Klaid^a) + w(^w)] + 

- ^(d^j a )((T*s v<7 + CJS va ) |. 

This expression contains no explicit spinor terms, as required, but we can improve it by elimi¬ 
nating the rank-2 tensors s pi/ and by using the Fierz identity 

a^ = (a 2 - a; 2 )" Ve^ - u8 pvpa )j p k a , (141) 

V" = (a 2 - w 2 )- 1 (we w ' 1 - o^ vpr, )j p k a , (142) 

where we define the partially antisymmetric object 

5 ^ P a = _ rp’rf'p), (143) 


We also require the derivatives of these identities, which are 

Jl , , ,2 




id *s — 

'-'fi ’Jya — 


{[2mj(d p u;) - (a 2 + cj 2 )(d p a)}e uape + [2cuj(d M cr) - (a~ + u; 2 )(d p uj)]S uape }j p k e 


+ 


(a 2 -co 2 ) 2 

{atvcrpe ~ u8 vapfL )[{d, l j p )k' L + j p (d M fc e )] 


(144) 


{[-2<7w(<9 p cr) + (cr 2 + w 2 )(9 Al w)]e i , CTpe + [-2erw(<9 p w) + (a 2 + u} 2 )(d p a)\S vcrpe }j p k e 


+ 


{a 2 - uj 2 ) 2 

(uje UC rpe ~ U S vape )[(d p j P )k e + j P {8 p k e )] 


(145) 


the rank-2 dependent terms in (114011 become 

-^j a W(dp*Sucr) + u{d„a va )] = ^(cr 2 - w 2 ) -1 {2 au}k„[u(d p cr) - a (<9 p w)] 

+ 2i cru}e u(rpe j rT (dp,j p )k e + j v j a (dp.k^ia 2 + w 2 )} - j^(d p k„)(a 2 + co 2 ), (146) 

- ^{dp.f)(o*s va + us va ) = {a 2 - w 2 ) -1 ( j^jifj a (d p k a )(a 2 + w 2 ) - ^-aue vape (d p ,j ,T )j p k t 


12 ^ 


6 


- — fc^cr 2 + w 2 )[w(<9 p u;) - cr(<9 p cr)] 


(147) 


giving us the final form of our identity 

[^{dpip) - (5^)7^] = (cr 2 - w 2 ) _1 {fc y [w(5 p CT) - cr(a M w)] - it„ape{d p f)j p k e 
- i j v m a (dp,n a )}. 


(148) 
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